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The elsewhere surmised topological origin of phase transitions is given here new important evi- 
dence through the analytic study of an exactly solvable model for which both topology and thermo- 
dynamics are worked out. The model is a mean-field one with a fc-body interaction. It undergoes 
a second order phase transition for k — 2 and a first order one for fc > 2. This opens a completely 
ri , new perspective for the understanding of the deep origin of first and second order phase transi- 

O ' tions, respectively. In particular, a remarkable theoretical result consists of a new mathematical 

D ' characterization of first order transitions. Moreover, we show that a "reduced" configuration space 

can be defined in terms of collective variables, such that the correspondence between phase transi- 
tions and topology changes becomes one-to-one, for this model. Finally, an unusual relationship is 
t^ ' worked out between the microscopic description of a classical A'^-body system and its macroscopic 

4-J , thermodynamic behaviour. This consists of a functional dependence of thermodynamic entropy 

. ■ upon the Morse indexes of the critical points (saddles) of the constant energy hypersurfaces of the 

microscopic 2A'^-dimensional phase space. Thus phase space (and configuration space) topology is 
directly related to thermodynamics. 



PACS numbers: 05.70.Fh; 02.40.-k; 75.10.Hk 



I. INTRODUCTION 



Thermodynamical phase transitions are certainly one of the main topics of statistical physics. A huge amount 

of work, both theoretical and experimental, has been done during the past decades leading to remarkable successes 

r^ _ as is witnessed, for example, by the Renormalization Group theory of critical phenomena. However, there are still 

\^ . longstanding open problems about phase transitions; among them we can mention amorphous and disordered systems 

^^ ■ (like glasses and spin-glasses) undergoing "dynamical" transitions, or first-order phase transitions, which are still 

\l lacking a satisfactory theoretical understanding of their origin. Moreover, on the forefront of modern research in 

^^ statistical physics, standard theoretical definitions and methods are challenged by the experimentally observed phase 

transitions occurring in small classical and quantum systems (nano and mesoscopic) like atomic and molecular clusters, 

polymers and proteins, Bose-Einstein condensates, droplets of quantum liquids, etc. Finally, in the mathematically 

rigorous background of phase transitions theory, neither in the Yang^ee theory for the grandcanonical ensemble [l| 

nor in the Ruelle, Sinai, Pirogov theory for the canonical ensemble 0, an a-priori mathematical distinction can be 

(~| ' made among the potentials leading to first or second order phase transitions, respectively. 

O The present paper aims to contribute to the advancement of a recently proposed theoretical framework where the 

O singular behaviours of thermodynamic observables at a phase transition are attributed to major topology changes in 

J> ', phase space and - equivalently - in configuration space Qj S IE Ifl ■ More precisely, in Ref . @j S S S ' i* has 

r '~j • been proposed that thermodynamic phase transitions could be a consequence of suitable topology changes of certain 

rS ■ 
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submanifolds of the configuration space defined by the potential energy function. The presence of such topological 
changes has been recently shown to be a necessary condition, under fairly general assumptions, for the presence of 
a phase transition [SJ; however, the converse is not true, and no rigorous results are available yet on the sufficient 
conditions. The analytical or numerical study of particular models thus remains crucial to get hints towards more 
general results (see also Refs. ^ [l^ for recent results on one-dimensional systems, and Ref. [ll| as to the fully 
connected spherical model). 

In this perspective, we stress that the topological approach has been hitherto applied only to systems undergoing 
second-order phase transitions. Therefore tackling first order phase transitions in the topological framework is of 
great interest, because some new insight into the challenging problem of their origin can be obtained, and because 
this reinforces the working hypothesis that the topological approach could unify the treatment of the different kinds 
of phase transitions in view of encompassing also more "exotic" ones, like glassy transitions and the others mentioned 
above. 

In this paper we first study a model which, according to the value of a parameter, has no transition or undergoes a 
first or a second order transition. Remarkably, for this model an exact analytical computation is possible of the Euler 
characteristic - that is a topological invariant - of those submanifolds of the configuration space whose topological 
changes are expected to be related to the phase transition. We find J2J that the phase transition is actually signaled 
by a discontinuity in the Euler characteristic x{e) and that the sign of the second derivative of x(e) indicates the order 
of the transition. Then, we study the topology of submanifolds of a "reduced" configuration space, i.e., the space of 
some collective variables: in this case we find a one-to-one correspondence between phase transitions and topology 
changes. 

Finally, we derive a general result showing that an analytic estimate of another topological invariant of the same 
submanifolds can be worked out which allows to directly link topology and thermodynamic entropy. 

II. A KEY STUDY 

In this Section we present a study of the thermodynamical properties of the mean-field fc-Trigonometric Model 
(fcTM), as well as of the topological properties of its configuration space. A preliminary study of this model along these 
lines has already been reported in 12] ; there only the microcanonical thermodynamics was considered, while here we 
are going to discuss also the canonical thermodynamical properties. Being a model with long-range interactions which 
may undergo also first-order phase transitions, we expect canonical and microcanonical thermodynamic functions to 
be different, at least close to first-order transitions 



The fcTM is defined by the Hamiltonian: 



^ 1 
=1 



Hk=}_^^7Tf+Vkicpi,...,ipN) (1) 



where {ipi} are angular variables: (pi G [0, 2tt), {iTi} are the conjugated momenta, and the potential energy V is given 

by 

^k = j;^j-T Y. [l-cos((/.,, +... + <^,J] , (2) 

ii,...,ifc 

where A is the coupling constant. In what follows only the potential energy part will be considered. This interaction 
energy is apparently of a mean-field nature, in that each degree of freedom interacts with all the others; moreover, 
the interactions are fc-body ones. 

The fcTM is a generalization of the Trigonometric Model (TM) introduced by Madan and Keyes [Ij| as a simple 
model for the Potential Energy Surface (PES) - the hypersurface defined by the potential energy as a function of the 
N degrees of freedom - of simple liquids. The TM is a model for N independent degrees of freedom with potential 
energy ^ with k = 1: Vk=i. 

It shares with Lennard-Jones like systems |15| the existence of a regular organization of the critical points of the 
potential energy above a given minimum (the elevation in energy of the critical points is proportional to their index) 
and a regular distribution of the minima in the configuration space (nearest-neighbor minima lie at a well defined 
Euclidean distance). The PES of the kTM maintains the main features of the TM |lg, introducing however a more 
realistic feature, namely the interaction among the degrees of freedom (in the form of a /c-body interaction). 

Using the relation 

cos((^,, + ... + if, J = i?e(e''^-i ■ . . . ■ e^'^'O , (3) 



the configurational part of the Hamiltonian can be written as 

Vfc = NA[l-Re{c + is)''] 

[k/2] 



NA 



n=0 



(4) 



where c and s are collective variables, functions of {(pi}: 



C = —^COSipi , 
i 

s = j^'^sinip, . (5) 



We observe also that the model has a symmetry group obtained by the transformations 



ip^^ip^+ i^ 

ipi ^ -ipi (6) 



If we think of tpi as the angle between a unitary vector in a plane and the horizontal axis of this plane, we find that 
the first transformations are rotations in this plane of an angle i^ and the second is the reflection with respect to 
the horizontal axis. This group is also called Ckv 

Let us now derive the thermodynamical properties of the fcTM. 

A. Canonical thermodynamics 

The partition function is 

Zk = fd{ip}e-^"^ = |d{^}e~'5^^[i-«^(-+-)'l ; (7) 

introducing (5- functions for the variables c and s, 

and using the integral representation of the (5-function, we obtain for Zk 

Zk = fd{ip} fdxdy f N^^^ ^^NX(.-c) ^^N^.(y-s) g-/3AJV[l-fle(.+.y)'=I ^ 

j dx dy e-/5^^[l-«-(-+»I')1 f N^^^ e^N{^x+^.y) j ^y^^ ^-.NiXc+^s) _ (9) 

The last integral is easily computable using Eq. ((SJ, 

f d{ip} e-'^^^''+>"'^ = ( f d^ e-*(^™"'^+''"'"'^M , (10) 

and can be written in term of the Bcssel function Jq 

r2TT 

dip e-^Acos(>^-v) = 2tt Jo(A) , (11) 





where A = A cosip, fi = A smip, and A = v A^ + /i^. The partition function can then be written as 

Zk = N^{2nf~^ [ dx dy d\ d^i e-^'(-»^-»A<y+'3A-/3Ai^e(x+,:y)^-iog(Jo(A)) ^ ^^3) 



and since the Jq function is always positive, there are no problems in defining its logarithm. 

We want now to perform a saddle-point evaluation of the integral, so we have to look for the minima of the exponent 
in the complex A, fi, x, y plane. We note that if that points do not lie on the imaginary axis of the A, /i planes, the 
free energy of the model would be imaginary. So wc can safely rotate the integration path on the imaginary axis in 
the A, /i planes, which corresponds to the substitutions: zA — > A and i^ — *■ fi, then lA — > A and Jo (A) — > /o(-'^)j where 
Jo is the modified Bessel-function: 



''''''T.I 



27T 

A cos if 



dip e^™^'^ . (13) 



Zk = N\27t)^-^ / dx dy dX dii e-^9'=(^.2'.^.'^;/3) ^ (^4^ 



In conclusion we obtain 



where gk is the real function 

gk{x, y, A, m; /?) = /3A - Aa; - /xy - /3A Re{x + lyf - log(/o(A)) . (15) 

In order to find the stationary points, we first determine the subspace defined by the equations 



|i = 0, ,17, 



obtaining the relations 



A = -~(i /:^kRe{x + iyf-^ , (18) 

fi = f3 A k Im{x + iy)''-\ (19) 

thus we get 

A:^(3Ak\{x + iyf-^\. (20) 

Now, using Eqs. (|18I19|I . we can substitute A and /i with x and y in Eq. H15(l . obtaining, in term of the complex 
number z = x + iy, 

5fc(z;/3) = /3A + /3A(fc-l)i?ez'=-log/o(/3Ap|fcP-i|) , (21) 

and using the polar representation z = pe^^ 

gk{p,t,P) = l3A + l3A{k-l)p'' cos{kij)-\ogIo{l3Akp''-^). (22) 

The derivative with respect to ip leads to 

-/3A(fc - l)fcp'= sin(/fcV) = , (23) 

so that there are 2fc solutions 

TiTT 

V'n = -^ (n=l,...,2fc). (24) 

Observing that cos(A:'0„) = (—1)" we obtain 

gk{p,n; /?) = /? A + (-l)"/3 A (A: - 1) / - log/o(/3Afcp'=-i) (25) 

and we can restrict ourselves to n = 0, 1. Finally, the derivative with respect to p leads to the stationary points 
equation: 



hiPAkp''-^) 
loiPAkp"-^) 






where the modified Bessel-function Ji is defined by 

/i(A) ^^ r d^ cos^ e^™^'^ = /^(A) . (27) 

For n = 1 we have only the trivial solution p ~ 0, because the / functions are always positive. By using an expansion 
for small p one can show that this solution is a maximum for g. So we can study only the case n = 0. We note that 
if there is a non trivial solution (i.e., p{(3) / 0) of Eq. H26() . calling gk{(3) the value of gk{(3, /5(/3)), we have 

Zk ^ N^ {27r)''-^ e~^^U0K (28) 

and the free energy and internal energy are, respectively, 

fk{f3) = /3-'5fc(/3)-r'log(27r), (29) 

eu{P) = A(l - p^) . (30) 

Let us now analyze the case k = 1. In this case the solutions p — Q are not present, so that we have only the 
solution 

There is no phase transition, and using Eq. (|30|l we have 

This is the free energy of trigonometric model that has been mentioned before. 

For fc = 2 the solution p = is stable for high temperatures, but a non trivial solution of Eq. H2()|l appears at 
/3A = 1. The transition temperature is given by the condition 



<P92{p\Pc) 



dp2 



= 2/3cA(l - (3cA) = , (33) 

p=0 



so that we obtain f3cA. — 1; the transition is continuous, and the order parameter is p. It is easy to show that x — (c) 
and y — (s) (e.g., by adding an external field of the form —N{hc + ks) to the Hamiltonian and performing the limit 
ft, fc — > 0); then the vector (i,y) is the mean magnetization of the spins represented by the ipi. As p is the modulus 
of the magnetization, for /3A > 1, when p 7^ 0, the C21, symmetry is broken. 

When fc > 2, the non trivial solution of Eq. (|25|l appears at a given (3' but becomes stable only at f3" > /3', so that 
/5(/3) and e(/3) are discontinuous at /3"; instead of the instability region fi' < f3 < /?", in the microcanonical ensemble 
a region where the specific heat is negative appears, as we shall see below. The Ckv symmetry is broken in the low 
temperature phase, so that p can be used as an order parameter in revealing the symmetry breaking, even if it is not 
continuous at (3" . The transition is then of first order, but keeps the symmetry structure of a second order one, i.e., 
in the low temperature phase there are fc pure states related by the symmetry group also in the case of the first order 
transition. 

In Fig. n we report the caloric curve, i.e., the temperature T = (3^^ as a function of the average energy (per degree 
of freedom) e, for three values of fc, fc = 1, 2 and 3. As previously discussed, the temperature is an analytic function 
of e for fc = 1; for fc = 2 the system undergoes a second order phase transition at a critical temperature Tc = A, that 
changes to first order for fc > 2. 

In Figs. 121 and 121 we report the average potential energy w as a function of the average energy e and the temperature 
T as a function of v, respectively. It is apparent that, for fc > 2, the phase transition point always corresponds to 

Another feature which shows up in Figs. |21and|31is that the average potential energy v never exceeds the value 
A, i.e., although the maximum of V/N is equal to 2A, the region w > A is not thermodynanucally accessible to 
the system. The reason for this is in the mean-field nature of the system and in the fact that we are working in 
the thermodynamic limit N — > cx). According to Eqs. Q, the potential energy can be written as a function of the 
collective variables c and s defined in Eqs. ((SJ, which are the components of the function whose statistical average is 
the order parameter, i.e., the "magnetization". In the thermodynamic limit these functions become constants, whose 
value coincides with their statistical average, and since (c) = (s) = for T > Tc, and from Eqs. I@J this implies v ~ A 
for all T>Tc. 

As we shall see below, this fact remains true also in the microcanonical ensemble, which, however, is not equivalent 
to the canonical ensemble for the present model, due to the long-range nature of the interactions. 
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FIG. 1: Temperature T as a function of canonical average energy e for three different values of k; for fe=l there is no phase 
transition, while for k—2 there is a second order transition and for fc > 2 a first order one. 



B. Microcanonical thermodynamics 



As in other simple mean field models, also in the case of the fcTM it is possibile to perform a calculation of the 
microcanonical partition function, or microcanonical density of states in phase space, given by 



^N.kiE) 






S{Hk ~ E) 



(34) 



The computation of fi is similar to that of Z in the canonical case, so that we will go through it with less detail. 
Using the integral representation of the delta function, we get 



^nAE) 



dp f d nd 'P ,0(H,^E) 
27r / iV! 



(35) 



Now, as we are looking for a saddle-point evaluation of the integral over /3, we can rotate the integration path on the 
imaginary axis in the complex-/? plane. This is justified because, as in the canonical case, the saddle-point is located 
on this axis. We can now perform the integration over the momenta and use the fact that Vk{(p) = Vk{c{ip), s{ip)), see 
Eq. Q), to obtain 



QnAE) = Cn p'^ J dp d^ dr, /3-f e^^^-^''^'")) • 
■Jd^^d[Ni^~c{ip))]S[N{r^~s{^))], 



(36) 



where p — N/L and the constant Cn gives only a constant contribution to the entropy per particle, i. e., it is at most 
of order e^. The last integral can be evaluated using again the integral representation of the delta function, and 
rotating then the integration path as previously discussed; it turns out to be: 



dp. dv 



-N{pe,+vi-i) 



d"^e 



^ ^ (/J. COS t^i+z-' sin i^i) 
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FIG. 2; Canonical average potential energy w as a function of canonical average energy e for fc=l, 2 and 3. The upper phase 
transition point is, for Vfc > 2, Uc = A. 



^e'^(^«+-')(27r/o(A))^ , 



having defined A = yj [J? + v^-^ /q is a Bessel function as before. We can then write the density of states as 

where m = (/3, ^, ri,fi,v), e — E/N and 

/fe(m, e) = /3e - /3A[1 - Rc(e + zt?)"] 



(37) 



Then, using the saddle-point theorem, the entropy per particle, s = S/N, is given by (fcs = l): 

Sfe(e) = lim — logilAr.fc(e) = max/fe(m, e) . 



(38) 



To find the maximum of /^(m, e) one can calculate analytically some derivatives of / to obtain a one-dimensional 
problem that can be easily solved numerically with standard methods. 

As already done in the case of the canonical ensemble, in Fig. 0] we report the microcanonical caloric curve, i.e., the 
temperature T as a function of the energy (per degree of freedom) e, T{e) = [ds/de]^-^ for three values oi k, k — 1, 
2 and 3. As in the canonical case, the temperature is an analytic function of e for fc = 1, while for fc=2 the system 
undergoes a second order phase transition at a certain energy value Cc, that changes to first order for fc > 2. 

We note that, for fc > 2, in a region of energies smaller than the critical energy Cc of the first-order phase transition 
the curve T{e) has a negative slope, i.e., the system has a negative specific heat. The fcTM is then another physical 
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FIG. 3: Temperature T as a function of canonical average potential energy v for three different values of k. 

model vifhere this feature is found (see, e.g., Ref. |lj| for other examples). This is not surprising at all since we are 
considering the microcanical thermodynamics of a system with long-range interactions; such a region is not present 
when we consider the canonical ensemble, as shown above; there, the region of negative specific heat corresponds to 
the region of instability of the non-trivial solution of the saddle-point equations. 

In Figs. 13 and we report the average microcanonical potential energy w as a function of e and the microcanonical 
temperature T as a function of v, respectively. It is apparent that, for k > 2, the phase transition point always 
corresponds to Wc = A. 

As in the canonical case, the average potential energy v never exceeds the value A, i.e., the region w > A is not 
thermodynamically accessible to the system also in the microcanonical ensemble. 



C. Topology of configuration space 

In this Section we want to investigate the relation between the phase transitions occurring in the /cTM and the 
topology of its configuration space. Given the potential energy V, the following submanifolds of configuration space 
are defined: 

My = {ip\ V{ip) < Nv} . 

As V varies from the minimum to the maximum allowed value of V{ip)/N , the manifolds M„ progressively cover the 
whole configuration space. These submanifolds, or their boundaries S^ = dMy, are then a possible way to depict 
the potential energy landscape of the system. The topology of the M^'s can be studied using Morse theory [l7l| : 
whenever a critical value of V{ip)/N - i.e., a value corresponding to one or more critical points, where the differential 
of V{(p)/N vanishes - is crossed, the topology of the Mv's change. It has been conjectured 3, 6] that some of these 
topology changes are the "deep" cause of the presence of a phase transition. The correspondence between major 
topology changes of the My 's and S„ 's and phase transitions has been checked in two particular models 0, Q ; more 
recently, it has been proved 8] that a topology change is a necessary condition for a phase transition under rather 
general assumptions ,28] . but the sufficiency conditions are still lacking. 
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FIG. 4: Microcanonical temperature T as a function of energy e for three different values of k; for k=l there is no phase 
transition, while for k=2 there is a second order transition and for fc > 2 a first order one. 



A natural way to characterize topology changes involves the computation of some topological invariants of the 
manifolds under investigation. One of such invariants is the Euler characteristic x'- ii^ ^; -2.] it was shown that the 
Euler characteristic x(v) of the submanifolds My and/or of the S^, shows a singularity in correspondence of the 
potential energy value Vc = v{ec) at which the transition takes place: this seems then to signal a particularly "strong" 
topology change. Remarkably, the Euler characteristic of My can be calculated analytically in our model. The general 
definition is flSl: 



N 



civ) = x(M.) = ^(-1)V.(^-) , 



(39) 



4 = 



where the Morse numbers ^i{My) are the number of critical points of index i of the function V{ip) that belong to 
the manifold My. As already mentioned, the critical points (called saddles in other contexts, e.g., in the physics of 
glasses) If are defined by the condition dVJt((^) = 0, and their index i (otherwise called the order of the saddle) is 
defined as the number of negative eigenvalues of the Hessian matrix 



nUv) 






(40) 



To be valid, Eq. H39|l requires that V{(p) is a Morse function, i.e., that its critical points are nondegenerate: this 
means that all the eigenvalues of the Hessian are nonzero at a critical point and that the critical points themselves 
are isolated. 

To determine the critical points we have then to solve the system 






= 



Vj-1,. 



,N 



that is, inserting Eq. |0J in the equations above, 

-A k Re[i{c + isf^^e'"^^] = A fc C''"^ sin[(fc - l)^' + (fij] = Vj 



,N. 



(41) 



(42) 
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FIG. 5: Microcanonical average potential energy v as a function of energy e for k=l, 2 and 3. The phase transition point is, 
for Vfc > 2, uc = A. 



where we defined c + is = Ce*'^. From Eq. Q we have 

Vk{ip)^NA[l-C''cos{ki^)] 



(43) 



then all the critical points with C(c^) = have energy v — V{(p)/N — A. We note that they correspond to vanishing 
magnetization. Let us now consider all the critical points with C,{(p) ^ 0. Then Eq. 142|) becomes 



and its solutions are 



sin[(fc - 1)V' + (^j] = Vj = l,...,7V. 



(pf = [mj-K - {k - l)V']mod 27r , 



(44) 



(45) 



where rUj G {0, 1}. Since in Eq. (|42|) ( appears to the k — 1-th power, in the case fc = 1 Eqs. (|42|) and H44|) coincide. 
This means that the solutions given in Eq. I|45|l are all the critical points, regardless of their energy, in the case k — 1 
and all the critical points but those with energy ?; = A in the case fc > 1. The critical point (^™ is then characterized 
by the set m = {nij}. To determine the unknown constant ip "we have to substitute Eq. H45|) in the self-consistency 
equation 



3 3 

If we introduce the quantity n{(p) defined by 



n = N ^ 2, "^i 



which means 



1 - 2n = iV-i 2J(-1)' 

3 



(46) 



(47) 



(48) 
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FIG. 6: Microcanonical temperature T as a function of microcanonical average potential energy energy v for three different 
values of k. 



we have from Eq. (|46|l 



C=|l-2n| 



tpi 



2lTr/k 

{2l + l)Tr/k 



for n < 1/2 , 

for n > 1/2 , 



(49) 
(50) 



where I £ Z. Then the choice of the set {rrij} is not sufficient to specify the set {(fij}, because the constant ip can 
assume some different values. This fact is connected with the symmetry structure of the potential energy surface: the 
different values of V'z correspond to the symmetry-related critical points under the group Ckv ■ 
We can then state that all the critical points with ( ^ - whose energy v ^ A - have the form 



The Hessian matrix is given by 






(pf' = [rrijTT - (fc - l)V'i]mod 2Tr ■ 

A k Re[N-^{k - l)(c + is)fc-2e»(v,+¥'j) 



In the thermodynamic limit it becomes diagonal 

Sij A k (''^^ cos (■0(fc - 1) + ifi) 






(51) 



(52) 



(53) 



One can not a priori neglect the contribution of the off-diagonal terms to the eigenvalues of 7i, but we have numerically 
checked that their contribution change at most the sign of only one eigenvalue over N: we note that in the case of 
the mean- field XY model this result has been proven explicitly J>|- Neglecting the off-diagonal contributions, the 
eigenvalues of the Hessian calculated in the critical point (f are obtained substituting Eq. 1)5 l|l in Eq. ()53|1 : 



Xj = (-1)"'^A k C'^-i 



(54) 
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so the index of the critical point is simply the number of rrij = 1 in the set m; we can identify the quantity n{ip) given 
by Eq. H47|) with the fractional index v/N of the critical point (p. Then, from Eq. Q, H49|l and H5U|) we get a relation 
between the fractional index n{if>) and the potential energy v{lp) = V{ip)/N at each critical point ip: 



n{v) 



1 




/_, 


V \ 




V 


i/fc- 





1 - 


- sgn (^1 - 




1 - 


- 




2 




" A/ 




A 





(55) 



Moreover, the number of critical points of given index i^ is simply the number of way in which one can choose v times 
1 among the {"nij}, see Eq. (|51ll . multiplied for a constant Ak that takes into account the degeneracy introduced by 
Eq. |5n|). 

We have then completely characterized the critical points with ( =/= 0. Now we are going to argue that, in order 
to compute the Euler characteristic of the manifolds Aly, we can neglect the critical points with C = Oi thus showing 
that the knowledge of the critical points considered so far is sufficient. The critical points with C — a-''^ degenerate: 
the Hessian vanishes at these points. This means that the potential energy is no longer a proper Morse function when 
V > A, and we could use its critical points to compute the Euler characteristic of the manifolds My only when ?; < A. 
To overcome this difficulty we reason as follows. Morse functions are dense in the space of smooth functions: this 
means, in practice, that if a function / is not a Morse function, any small perturbation will make it a proper Morse 
function [19j, and we can consider, as our Morse function, the function Vk, i.e., the potential energy plus a linear term 
which can made as small as we want: 



N 



Vk = Vk + Y, h. 



'-Pi 



(56) 



where h £ M^. The perturbation changes only slightly the critical points with C 7^ 0) but completely removes the 
points with C = for any h y^ 0, no matter how small. All the critical points of this function are given by the solutions 
of the equations 



sin[(fc - 1)^ + <fj] = hj V j = 1, . . . , iV , 



(57) 



which are only a slight deformation of Eqs. (|44|l . so that provided all the /I's are very small the numerical values of 
critical points and critical levels will essentially coincide with those computed so far, in the case /i = but assuming 

The fractional index n = v/N of the critical points is a well defined monotonic function of their potential energy w, 
given by Eq. (|55|) . and the number of critical points of a given index v is ^fc(^)- Then the Morse indexes ^^{My) of 
the manifold My are given by Aki^,) if v/N < n{v) and otherwise, and the Euler characteristic is 



Nn{v) 

x{v)=Au E(-l)' 

v=0 



-^^(-)^"'"^ (!(.;)' 



(58) 



using the relation Tll^i-^^ O = (-1)"(^™')- 

In Fig.[7|we plot cr(w) = hmjv^oc -^ log |x(w)|, that, from Eq. (|^ . is given by: 

a{v) — ~n{v) logn(u) — (1 — n{v)) log(l — n(v)) 



(59) 



It has to be stressed that cr(ti) is a purely topological quantity, being related only to the properties of the potential 
energy surface defined by 14 ((p), and, in particular, to the energy distribution of its saddle points. From Fig. [7|we 
can see that there is an evident signature of the phase transition in the analytic properties of cr{v). First, we observe 
that the region w > A is never reached by the system, as discussed before and showed in Figs. [S] and El as to the 
microcanonical case, and in Figs. El and El as to the canonical case; this region is characterized by a-'{v) < 0. The main 
observations are that: i) for fc=l, where there is no phase transition, the function a{v) is analytic; ii) for fc=2, when 
we observe a second order phase transition, the first derivative of cr(v) is discontinuous at Vc = v(ec) — A, and its 
second derivative is negative around the singular point. Hi) for fc > 3 the first derivative of <j{v) is also discontinuous 
at the transition point Vc — A, but its second derivative is positive around Vc- In this case a first order transition 
takes place. Therefore the investigation of the potential energy topology, via (y{v), allows us to establish not only the 
location but also the order of the phase transitions, without introducing any statistical measure. 

The previous results suggest us to conjecture that there is a relation between the thermodynamic entropy of the 
system and the topological properties of the potential energy landscape, as probed by ct(u). We recall that the 
presence of a first order transition with a discontinuity in the energy is generally related l2y| to a region of negative 
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FIG. 7: Logarithmic Euler characteristic of the M„ manifolds (^{v) (see text) as a function of the potential energy v. The phase 
transition is signaled as a singularity of the first derivative sd, Vc = A; the sign of the second derivative around the singular 
point allows to predict the order of the transition. The region ii > A, in which cr'{v) < 0, in not reached by the system (see 
text). 



specific heat, i.e., of positive second derivative of the entropy. Thus, it seems that at least around the transition point 
the thermodynamic entropy and a{v{e)) are closely related: more precisely, it seems that the jump in the second 
derivative of s(e) is determined by the jump in the second derivative of cr(w(e)). Then it should be possible to write 



s(e) -a(^;(e))+7^(e) 



(60) 



where TZ{e) is analytic (or, at least, C^) around the transition point. 

In Sec. lIIII we explain how such a hnk between thermodynamics and topology could be obtained. But before doing so 
we show a different way of looking at the topology of the submanifolds of configuration space defined by the potential 
energy. 



D. Topology of the order parameter space 



A feature of many mean-field models (although not of all of them) is that the potential energy can be written as a 
function of a collective variable, whose statistical average is the order parameter. In the case of the fc-trigonometric 
model this variable is the two-dimensional "magnetization" vector defined as m = {c,s), where (see Eqs.jS)) 



i=l 

1 ^ 



(61) 



i=l 
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FIG. 8: The submanifolds A4v in the case k — 1 for v — 0.5A, A, 1.5A, 2A (from left to right). All the submanifolds are 
topologically equivalent to a single disk. 
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FIG. 9: The submanifolds Aiv in the case k = 2 for v = 0.5A, A, 1.5A, 2A (from left to right). As v < Vc = A the submanifolds 
are topologically equivalent to two disconnected disks, while as v > Vc they are equivalent to a single disk. 

Written in terms of (c, s), the potential energy is a function defined on the unitary disk in the real plane, which is 
given by (see Eq. Ql) 



Vkic,s)^NA 



[k/2] 



1-E 



n=0 



2n 



(-1)" C 



n ^k — 2n 2n 



In the particular cases k = 1,2,3 the potential energy Vk reads as 

Vi{c,s) = NA{l-c) , 
V2{c,s) = 7VA(1 -c^ + s^) , 
V'3(c,s) = NA{l-c^ + 3cs^) 



(62) 



(63) 
(64) 
(65) 



and it is then natural to investigate the topology of the My^s seen as submanifolds of the unitary disk in the plane, 
i.e., we now consider the submanifolds 



M„ = {ic,s)eD^ I Vk{c,s)<Nv}. 



(66) 



where D^ = {{c,s) e R^ | c^ + s^ < 1}. The 7Wt,'s are nothing but the Af^'s projected onto the "magnetization" 
plane. 

The topology of these manifolds can be studied directly, by simply drawing them. In the case k — 1, where no phase 
transition is present, no topology changes occur in the 7Wt,'s, i.e., all of them are topologically equivalent to a single 
disk D^ (Fig.©- When k = 2,3, and a phase transition is present, there is a topology change precisely at Vc = A, 
where k disks merge into a single disk (see Figs. l^and llU|) . The detail of the transition, i.e., the number of disks 
which merge into one, clearly reflects the nature of the symmetry breaking for the particular value of k considered 
(similar pictures are obtained for fc > 3). Thus, when projected onto the order parameter space, the correspondence 
between topology changes and phase transitions becomes one-to-one (this was already found in 21] for the mean-field 
XY model); however, at variance with the study of the topology of the "full" M^s, no direct way to discriminate 
between first- and second-order transitions seems available in this picture. 



III. TOPOLOGY AND THERMODYNAMICS: A DIRECT LINK 



We now show how it is possible to establish a general relationship between topolo gy a nd thermodynamics. This 
can be achieved by improving some preliminary results on the subject reported in |22ll23| . 
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FIG. 10: The submanifolds A^„ in the case fc = 3 for v = 0.5A, A, 1.5A, 2A (from left to right). As v < Vc = A the submanifolds 
are topologically equivalent to three disconnected disks, while as v > Vc they are equivalent to a single disk. 



Consider a generic classical system described by a Hamiltonian 

N 

2 



1 ^ 
H= -Y.pl + V{q) 



where q — {qi, . . . ^qn) and the symbols have standard meaning. Then consider the microcanonical entropy S{E) 
defined as 



SiE)='^\ogn,{E), 



where v = 2N — 1 , where we shall put fc^ = 1 and where 



n,{E) 



1 

m 



da 



NHW 



(67) 



(68) 



with l|ViJ|| — {J2iPi + [^i^(9)]^}^^^j i-6-, ^jy is the microcanonical density of states. Here S^ is the constant- 
energy hypersurface in the 2iV-dimensional phase space F corresponding to the total energy E, that is T,e = 
{{pi,...,PN,qi,---,qN) G T\H{p,q) = E}. 
From the general derivation formula |24| 



^l^ada^l^AHa)da , 



(69) 



where a is an integrable function and A is the operator 

V 



A{a) 



NH\\ 



\7H 

Wm 



the following result is worked out 



dn^{E) 

dE 



m 



\VH\\ \\VH\\ \N 



(70) 



where M* = V(ViJ/||ViJ||). M^ is directly proportional to the mean curvature Mi of Eb seen as a submanifold 
of M^ [23 according to the simple relation Mi = —Mi/(2N — 1). By integrating equation Eq. H7U|) we obtain the 
following equivalent expression for fl^{E) 



n,{E) 



Nl 



drj 



da M^ 



1 



. ||vi?|| ||vi7|| m 



dp, 



Me 



M^ 

Wm\ 



(71) 



and then, at large N, considering that the volume measure dfi concentrates on the boundary S^, we write 



^! Jm. '^^ llVff II 



{5E) r da M{ 

^vT h^ Wm Wm 
m 

m 



\wH\n f 



da 



(72) 
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where, in the last approximate replacement, we have used that ||ViJ|| is positive and only very weakly varying at 
large N. 

By means of Holder's inequality for integrals we get 









da 



the sign of equality being better approached when Af* is everywhere positive. Hence, using Eqs. (|68|l . (|71|l and H72|l 



n^iE) < [Q.iE)]^ ( — 



1 



da 



iVffl 



m 



*\N 



SE 



iVi/ll 



(73) 



and introducing a suitable deformation factor d(E) we can write 



[n.iE)]^ = 



d{E)6E 

(WVHW) 



Nl 



da 



NH 



\M- 



*\N 



(74) 



so that 



n,{E) 



2^[d{E)]^{5E) 
(llVi/||)^+i 



\N 



da \K, 



n{E)\ 



(75) 



A'l A.f*|A' - c^. J ^ L .-..^^f - N\\Kg\ +n{E), with ki, . . . ,KAr the principal curvatures 

of ^E, and Kg is the Gauss-Kronecker curvature of E^;, Kg = Y[i=i ^i- T^i^) is a remainder. Again we have used 
that ||V7J|| is only very weakly varying at larg e TV and that it is always positive. 
According to the Chern-Lashof theorem [2q 



where we have used 2 '^ \M^\'^ = (ki + K2 + • ■ • + kn) 



/ da\KG\^voliS^-^)y^^i,iEE) 



(76) 



where Sj^ ~ is an N — 1-dimensional hypersphere of unit radius and /ii(S]£;) are the Morse indexes of 'Se which are 
defined exactly as those of the M^'s seen in the previous Sections, i.e., as the number /i of critical points of index i on 
a given level set T,e = E[^^{E); a critical point is a point where Vi? = 0, the index i of a critical point is the number 
of negative eigenvalues of the Hessian of H computed at the critical point. 
Finally the entropy per degree of freedom reads as 



S{E) = ^logn^iE) 



N 



log 



N 

vol{S^-')Y,t^^i^E)+ da 

1=0 ■'^E 



2NN\ 



j_ 2^[d[E)]^{5E)^ 
N °^ (||VH||)^+i 



(77) 



The meaning of Eq. (|77|l is better understood if we consider that the Morse indexes ^i{M) of a differentiable manifold 
M are related to the Betti numbers bi{M) of the same manifold by the inequalities 



fi,{M) > h{M) 



(78) 



The Betti numbers are fundamental topological invariants |l8| of differentiable manifolds; they are the diffeomorphism- 
invariant dimensions of suitable vector spaces (the de Rham's cohomology spaces), thus they are integer numbers. 
The equality sign holds only for the so-called perfect Morse functions, which are rare, however, at large dimension we 
can safely replace Eq. H78|l with iii{M) ~ bi{M) (see e.g. Refs. 7]). 
Equation (|77|l . rewritten as 



S{E) ^ llog 



TV 

wo;(§f-i)V6,(I]£)+ / da 



4=0 



2^N\ 



^_ 2^[d{E)]^{5E) 
N °^ (||ViJ||)^+i 



(79) 
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links topological properties of the microscopic phase space with the macroscopic thermodynamic potential S{E). 

In particular, even though the function TZ{E) is unknown, sudden changes of the topology of the hypersurfaces S^ 
(reflected by the energy variation of Yl 6i(Es)) necessarily affect the energy variation of the entropy. 

Finally, we resort to the fact that - at large N - the volume measure of I^e concentrates on Sf^^!i/2 ^ ^^{y)j where 

^(2K^/2 = {(Pi' ■ ■ ■ 'Pn)\ J2p1 = C^K)} is the kinetic energy hypersphere and M(^v) = {(^i, ■ • ■,qN)\V{q) < (V)}, 
so that TiE can be approximated by this product manifold, and we resort to the Kunneth formula |18j for the Betti 
numbers of a product manifold A x B, i.e. 

h{A^B)= Y. bj{A)bk{B) (80) 

j+k=i 

which, applied to S^^u/2 x ^(y), gives ^^(Sb) = 2b,{My) for i = 1, . . . , TV - 1, and 6j(Eb) = bj{My) for j = 0, TV, 
since all the Betti numbers of an hypersphere vanish but bo and b]sf which are equal to 1 18]. Eventually we obtain 



S{v) c^ llog 



vol{§ 



i~')(ba+J2^b'{Mv) + bN] 

1 , 2^[d{E)]'^{SE)'^ 



+ niEiv))\ + - log -|| — 1^— , (81) 

where TZ{E{v)) stands for the integral on the product manifold of the remainder which appears in Ea. (|79|l . The 
equation above makes explicit the fact that the kinetic energy term of a standard Hamiltonian is topologically trivial. 

From this equation we see that a fundamental topological quantity, the sum of the Betti numbers of the submanifolds 
My = {((?!, . . . ,(7Ar) G M^|\^((7) < w} of configuration space, is related, although with some approximation, to the 
thermodynamic entropy of the system. 

While a relationship between topology and thermodynamics exists, as is shown by both Eqs. (|81|1 and H77|) . an 
analytic formula linking the Euler characteristic to thermodynamics has not been found yet and is unlikely to exist. 
Therefore, in those cases allowing the analytic computation of the Morse indexes (as for the models tackled in this 
paper) , besides the computation of the Euler characteristic through the formula H39|l , we can use the Morse indexes 

to replace the sum (&o + X]j=T 26i(Mt,) + &Ar j in Eq. (|^ with (/^o + X]j=T 2^i(Mt,) + //at j , having resorted again 
to the estimate bi{M) ~ fii{M). Then we can plot this sum, that we shall denote by /i, against the potential energy 
density for the three different cases considered: k — 1,2, 3. The result is shown in Fig. ^2 where sharp differences are 
evident between the three situations: i) k — 1, absence of any phase transition, in which case the pattern of /i vs v is 
smooth; ii) k = 2, second order phase transition, in which case the pattern oi pi vs v displays a sharp non-differentiable 
change at the phase transition point; Hi) k = 3, first order phase transition, in which case the pattern oi fi vs v again 
displays a sharp non-differentiable change at the phase transition point which is approached from the left with an 
opposite concavity with respect to the second order transition case. Likewise the Euler characteristic, the quantity fi 
is in one-to-one correspondence with topology changes of the manifolds My, but fi has the advantage of being directly 
related with thermodynamic entropy. 

Both the general analytic result of Eq. (|81ll and the particular analytic result obtained for the fcTM and reported 
in Fig. 111! are of great relevance in view of a deeper understanding of the relationship between topology changes of 
configuration space and phase transitions: further work is ongoing along this direction. 

As a final comment, let us remark that the clearcut differences among the three different cases in Fig. ^2 are 
detected prior to and independently of the definition of any statistical measure in configuration space: the relevant 
information about the macroscopic physical behavior is already contained in the microscopic interaction potential 
itself and concealed in its way of shaping configuration space submanifolds. 

IV. CONCLUDING REMARKS 

We have presented a mean-field model whose thermodynamics is exactly solvable in both the canonical and the mi- 
crocanonical ensemble: the model depend on a parameter k and exhibits no transitions, a continuous phase transition 
and a discountinuous one, in the cases k — 1, k — 2, and A: > 3, respectively. For this model, a clear and sharp relation 
between the presence of a phase transition and a particular topology change in the submanifolds of the A'^-dimensional 
configuration space is analytically worked out: this correspondence becomes one-to-one if we look at the submanifolds 
of a reduced two-dimensional configuration space spanned by collective variables. Moreover, not only the presence and 
the location in energy of the transition can be inferred by looking at the behavior of topological quantities: also the 
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FIG. 11: Logarithm of the sum of the Morse indexes divided hy N, fx = -^ log mo + 2^^^j^ fii{A4v) + fiN 
Mv versus the energy density v, scaled with A, for k = 1, 2, 3. 



of the manifolds 



order of the transition is related to the behavior of a topological invariant of the above mentioned submanifolds, their 
Euler characteristic. This suggests that topological quantities are linked in general to thermodynamic observables. 
Such a general link, although based on some approximations, has been derived in the final Section of the paper. 

The results presented here confirm the validity and the potential of the topological approach to phase transitions, 
which has recently received a rigorous background via the proof of a theorem |^ stating that, for systems with 
short-ranged interactions, topology changes in the submanifolds of configuration space are a necessary condition for 
a phase transition to take place. The model we studied here is not short-ranged, thus the theorem might probably be 
extended to a more general class of systems. However, we would like to mention the case of a mean-field model, the 
fully coupled ip^ model, which has been recently studied [23, and where the relation between the topology changes 
in the submanifolds of the configuration space and the phase transition is less straightforward. Further work is then 
needed to assess the full potential and the limits of the topological approach to phase transitions. 
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